Effect of a fermion on quantum phase transitions in bosonic systems 
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The effect of a fermion with angular momentum j on quantum phase transitions of a {s, d) bosonic system is investigated. It is shown 
that the presence of a fermion strongly modifies the critical value at which the transition occurs, and its nature, even for small and 
moderate values of the coupling constant. The analogy with a bosonic system in an external field is mentioned. Experimental 
evidence for precursors of quantum phase transitions in bosonic systems plus a fermion (odd-even nuclei) is presented. 
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Quantum phase transitions (QPT) are qualitative changes in 
the ground state properties of a physical system induced by 
a change in one or more parameters in the quantum Hamil- 
tonian describing the system. Originally introduced in the 
,1970's 111121, they have been the subject in recent years of many 
investigations and have found a variety of applications in many 
areas of physics and chemistry |i3i,|4i]. One of these applica- 
tions is to atomic nuclei, where ^Ts have been extensively 
investigated (for a review, see (E, 6,'?']) within the framework 
t)f the Interacting Boson Model (IBM), a model of even-even 
nuclei in terms of correlated pairs of valence nucleons with 
angular momentum J - 0,2 treated as bosons (s,d) |8]. For 
this case also finite size effects 1^ IT ] and scaling behav- 
ior 112, Tsi [14I1 have been investigated, both analytically and 
numerically, showing that precursors of QPT can be seen even 
for relatively small values of A^. QPTs have also been extended 
to excited states quantum phase transitions, that is qualitative 
changes in the properties of the system as a function of the ex- 
citation energy tl5il . In this letter we present results of an inves- 
tigation of the effect of a fermion on QPTs in bosonic systems. 
We do this in atomic nuclei by making use of the Interacting 
Boson Fermion Model (IBFM), a model of odd-even nuclei in 
terms of correlated pairs with angular momentum J - Q,2{s,d 
bosons) and unpaired particles with angular momentum J - j 
a fermions) I1I61I . As an illustration we take j - 11/2. We 
note, however, that our method of analysis can also be used 
for systems with other values of the fermion, j, and boson, 
J, angular momenta, for example the spin-boson systems dis- 
cussed in UtIi . the simplest case of which is a fermion with 
j = 1 12 {i.e., a single spin) in a bath of harmonic oscillator 
one-dimensional bosons of interest in dissipation and light phe- 
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nomena. QPTs in IBFM for selected orbits have been investi- 
gated by Alonso et al. |18, T^. Here we focus on the effect of 
a fermionic impurity on QPTs in bosonic systems. Our main 
results are that, (1) the presence of a single fermion greatly 
influences the location and nature of the phase transition, the 
fermion acting either as a catalyst or a retarder of the QPT, and 
(2) there is experimental evidence for quantum phase transitions 
in odd-even nuclei (bosonic systems plus a single fermion). 

Within the context of the geometric collective model of nu- 
clei, the effect of an odd particle on collective properties was 
investigated years ago in core-particle models. However, both 
our results are novel, since (i) in the deformed phase the effect 
of the fermion is of order and thus vanishes in the geo- 
metric limit — > 00. However, we explicitly show that the 
effect is large in transitional nuclei even for small and moderate 
values of the coupling strength of the fermion to the bosons; 
(ii) the experimental evidence for QPTs in odd-nuclei has not 
been presented earlier and we show one of the key signatures 
of QPTs in nuclei, the two-neutron separation energies. This 
quantity is discontinuous for a first order phase transition at the 
transition point (or it has a sudden jump for a finite system). As 
shown below, this jump is observed experimentally. 

To prove our point, we consider the Hamiltonian of a system 
of {s, d) bosons coupled (with a quadrupole interaction) to a 
single fermion with angular momentum j Iil6ll 



H^Hb + Hf + V, 
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Here = d' ■ d is the number operator for d-bosons, = 
(d'''xs+s''xdy^^+X(d'^xdy^^ andq - (a\xajf^\ are quadrupole 
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operators of bosons and fermion respectively, eo is the scale of 
the boson energy, ej is the energy of the single fermion and 
r the strength of the quadrupole Bose-Fermi interaction. The 
dot and cross indicate scalar and tensor products and the ad- 
joint operators for bosons and fermions are = (-)'' d-^ and 



(-y 



QPTs of the purely bosonic part of the 



Hamiltonian Hb have been extensively investigated 
There are two control parameters ^ and x- For fixed as one 
varies f , < ^ < 1, the bosonic system undergoes a QPT. The 
phase transition is first order for;^' + and becomes second or- 
der at = 0. No phase transition occurs as a function of x- In 
this article, we take;^' - which case the two "phases" 

of the system have U(5) symmetry = 0) and SU(3) symmetry 

= 1) iSl- The critical point, separating the spherical [U(5)] 
and axially-deformed [SU(3)] phases, occurs at — 1/2. 

A complete study of the properties of quantum phase transi- 
tions necessitates both a classical and a quantal analysis, and a 
consideration of other couplin gs o f fermions to bosons in ad- 
dition to quadrupole coupling . In order to emphasize the 
main features of the results, we report here only the classical 
analysis. This amounts to constructing the combined Bose- 
Fermi potential energy surface (Landau potential) and minimiz- 
ing it with respect to the classical variables. To this end, we 
introduce a boson condensate llst] 



\N-li,y) ^ -^\bl{fi,y)\' \0) , 



(3) 



bl = (1 -HyS2)-i/2[^ cos yc/J-HyS sin r(t/2+t/!2)/V2-H.s^], in terms 
of the classical variables /?, y. The expectation value of Hg of 
Eq. (|2| in the condensate is lllOll 
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We then evaluate the expectation value of Hf and Vbf in the 
condensate thus obtaining a fermion Hamiltonian 
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The matrix gm^mjiN', /?, y) is a real, symmetric matrix 
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where C^'" denotes a Clebsch Gordan coefficient. The 

j,mi:j,-m2 

eigenvalues e, and eigenvectors t/f, of the matrix g are the single- 
particle energies and wave functions of the fermion in the de- 
formed (j6, y) field generated by the bosons. For y = 0° (field 
with axial symmetry), x - ^ ^nd /3 small (i.e., neglecting ff) 
they were obtained years ago by Nilsson |23]. Fory 0°;;^^ = 
and j6 small, they were investigated by Meyer-ter- Vehn [24] . We 
have solved the problem in its generality and details are given 
in \2^. Here, we consider, for simplicity, the case y — 0° for 
which the eigenvalues are given in explicit analytic form ll25[l 
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where Pj = [(2; - l)j(2j + 1)(; + 1)(2; + 3)]-"^ The quan- 
tum number K = j, j-2, j has the physical meaning of 
the projection of the angular momentum on the intrinsic z axis 
of the condensate. 

Once the eigenvalues have been obtained, one can calculate 
the total energy functional (Landau potential for the combined 
Bose-Fermi system) 



Ei{N;/3,r,^,x;n 



EB(N-/3,r,^,x) 

+ Sj + ei{N-l3,y-x-S) . (8) 



This expression is the algebraic analog of the total potential en- 
ergy surface, obtained in the macroscopic-microscopic Strutin- 
sky procedure ll26ll . Minimization of Ei with respect to /? and y 
gives the equilibrium values fSg, ye (the classical order parame- 
ters) for each state. In the simple case of y = 0°, ;t' = ~ 
the total energy functional becomes 

V7 



(9) 



EK{N-/3-^-ry) = Eb\N-/5,Q-^,- — 

V7 

+ AK(N;P;- — ;r) 
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Minimization (-^ = 0) gives the equilibrium yS^ values shown 
in Fig. [1] as a function of the control parameter ^ of the bosonic 
phase transition. 

By comparing the top part of this figure (purely bosonic 
system) with the bottom part, one can see that the effect of 
the fermionic impurity is to wash out the phase transition for 
states with K - 1/2, 3/2, 5/2 and to enhance it for states with 

= 7/2, 9/2, 11 12. In other words, the fermion acts as a cat- 
alyst for some states and as a retarder for others. Also, when 
the coupling strength becomes very large, the minima for some 
large K, in the figure K - 11/2, shift to negative values (oblate 
deformation). In addition to this result, also known qualita- 
tively from particle-core models, the effect of the fermion is to 
move the location of the critical point even for small and mod- 
erate values of F. Physical values of F in the eiPm, esEu, 65Tb 



nuclei, where the phase transition occurs, are F = -0.125 II27I1 . 

It is interesting to note that the effect of the fermionic im- 
purity is similar to the effect of an external field with linear 
coupling on a thermodynamic phase transition investigated by 
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Figure 1: (color online). Equilibrium values of the total boson plus ferraion 
energy (classical order parameter) as a function of the control parameter f of 
the bosonic phase transition for various values of the coupling constant T in 
units of £o and N = 10. The curves for F correspond, from left to right, to 
states with K = 1/2, 3/2, 5/2, 7/2, 9/2, 1 1/2. The value F = Ogives the purely 
bosonic case. 



Landau and Lifshitz years ago [28, p.456]. They considered the 
potential 

(^{rj) = <Do +Arf- + Brf + aT], (10) 

where a is the strength of the coupling to the external field and rj 
a classical variable (order parameter). (The bosonic part of this 
potential Aif + Bif has only a second order transition). After 
a projective transformation - if which does not change 
the nature of the phase transition and some rearrangement, the 
IBFM potential can be written for small 77, in the form 



EKiii) = hQ+Arf + Cif + Bif + aK\2ii + -l-jz 



(11) 



where ok is the strength of the coupling for each K value. By 
comparing the couplings in Eqs. (10) and (11), one can see that 
the IBFM Ek (??) is more general than 0(77) since it has a linear 
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Fi gure 2; (color online). The total energy at the equilibrium value, (£/^)min 
Eq. (9), as a function of the control parameter f for an intermediate value of the 
strength of the Bose-Fermi interaction, F = -0.25, with £0 = 1 and N = 10. 
The dashed curve, labeled by Eg, is the corresponding energy of the purely 
bosonic system. 



and quadratic coupling, but for small rj (note that \ri\ < 1), the 
quadratic term is negligible and the two expressions become 
identical. Also, in the IBM the bosonic part has a cubic term 
leading to the possibility of a first order transition. The analogy 
with bosonic systems in an external field also suggests that our 
results apply to the study of phase transition in superconductors 
in the presence of magnetic fields. 

Finally, having computed the equilibrium values, one can 
compute the total energies £, (A^; Pi,,7e',^,x) which for the spe- 
cial case discussed here are given by Ek (N; j5e', ^; F) and shown 
in Fig. 121 One can see again by comparing Ek with the energy 
of the purely bosonic system, Eb, that there is an effect espe- 
cially close to the critical value, ^c- The effect is not so much in 
the total energy Ek (where it is of order 1 jN) but in the deriva- 
tive of the total energy with respect to the control parameter, 

dEK 
flf ■ 

An important property of atomic nuclei is that they provide 
experimental evidence for shape QPTs, in particular, of the 
spherical to axially-deformed transition (U(5)-SU(3) symme- 
try) LSfct^J]. Three signatures have been used to experimen- 
tally verify the occurrence of shape phase transitions in nuclei: 
(a) the behavior of the order parameter (fie) as a function of the 
control parameter, measured through the B(E2) values propor- 
tional to PI; (b) the behavior of the ground state energies, mea- 
sured through the two-neutron separation energies, 5 2,1; and 
(c) the behavior of the gap between the ground state and the 
first excited 0^ state. Here for conciseness we concentrate only 
on S2n - -[Eq{N -I- 1) - £0(^^)1, which can be related to the 
derivative of the ground state energy, Eq, with respect to the 
control parameter, . S 2n can be written as a smooth contri- 
bution linear in the boson number A^, plus the contribution of 
the deformation Ulii] 
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In order to emphasize the occurrence of the phase transition 
it is convenient to plot the deformation contribution only, ob- 
tained from the data by subtracting the linear dependence, as 
a function of A^. In previous studies of the purely bosonic 
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Figure 3: (color online). The contribution of deformation to the two-neutron 
separation energies, S(2n)^cf for even-even 6oNd-62Sm-64Gd nuclei (top) and 
odd-even 6iPm-63Eu-(;5Tb nuclei (bottom), plotted as a function of neutron 
number. The contribution is enhanced in odd-even nuclei by approximately 300 
keV (at neutron number 92). Also the rise between neutron numbers 88 and 90 
is sharper in odd-even nuclei than in even-even nuclei. In the limit W — » oo 
(no finite size scaling) the quantity 5 (2«)def should be zero before the critical 



value and finite and large after that. The expected behaviour of - 
U(5)-SU(3) transition and W = 10 is shown in the inset. 



for the 



part it has been shown that is approximately proportional 
to the control parameter ^ fl. The experimental values of 
5(2n)def are shown in the top part of Fig. [3] for even-even 
nuclei (purely bosonic) and in the bottom part for odd-even 
nuclei (bosonic plus one fermion). They are obtained from 
the data fM] with A2,, = -14.61, -15.82, -16.997 MeV for 
Nd-Sm-Gd, respectively, and B2„ - 0.657 MeV, and with 
A2„ = -15.185, -16.37, -17.672 MeV for Pm-Eu-Tb, and 
B2n - 0.670 MeV. Precursors of the phase transition are visible 
in all six nuclei between neutron numbers 88 and 90 in both, 
and, most importantly, appears to be enhanced in odd-even nu- 
clei relative to the even-even case. 

In conclusion, we have presented here a classical analysis 
of quantum phase transitions in a system of bosons and one 
fermion (spin-boson system) and shown that (i) the addition of 
a fermion greatly modifies the critical value at which the phase 
transition occurs, and in some cases its nature; (ii) the effect 
is similar to that of adding an external field; (iii) there is ex- 
perimental evidence for these phase transitions in odd-even nu- 
clei at neutron number 88-90. The effect of the odd fermion 
is about 20% in 5(2n)def- A quantal analysis, in which the 



Hamiltonian H is diagonalized numerically for finite A^, pro- 
duces results similar to those of the classical analysis ll22ll . Our 
results are of interest not only for applications to nuclei, but 
also for applications to other systems in which a fermion is im- 
mersed in a bath of bosons, for example, the simple case of a 
spin 1/2 particle in a bath of harmonic oscillator bosons lll7n . 
Our analysis opens the way for a systematic study of QPTs in 
Bose-Fermi systems, in particular, of shape phase transitions in 
odd-even nuclei. This includes experimental studies and micro- 
scopic investigations using Density Functional Theory and/or 
other methods, in a way similar to what it has been done re- 
cently for the study of QPTs in purely boson systems (even- 
even nuclei) |,6,, 31J. 
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